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Abstract 



The standard classical expression of the state-resolved photo-dissociation cross section is not consis- 



tent with an efficient Bohr quantization of product internal motions. A new and strictly equivalent 
expression not suffering from this drawback is proposed. This expression opens the way to more 



realistic classical simulations of direct polyatomic photo-dissociations in the quantum regime where 

43. 
O , only a few states are available to the products. 
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I. INTRODUCTION 

Photo-dissociations play a key role in the evolution of planetary atmospheres (l| and inter- 
stellar clouds |2] and intense experimental research on their dynamics is currently performed 



0. 



4|. Their accurate theoretical description is thus a major issue 

The state-resolved absorption cross section measured in molecular beam experiments 
provides among the most detailed information on photo-dissociation dynamics [51] and we 
focus our attention on this quantity in the present work. 



Assuming that the electronic problem has been solved 



BE 



7|, state-of-the-art descriptions 



e performed within the framework of exact quan- 



of the previous observable are in princi 

turn treatments of nuclear motions [5|, l8l4l l| . However, despite the impressive progress of 
computer performances achieved in the last three decades, these approaches can hardly be 
applied to polyatomic processes as the basis sizes necessary for converging the calculations 
are usually prohibitive. 

An alternative approach is the quasi-classical trajectory (QCT) method [12J, |l3|, poten- 
tially applicable to very large molecular systems. Nevertheless, this approach obviously 
misses the fact that internal motions of the final fragments are quantized. This is a minor 
drawback in the classical regime where the number of states available to the products is 
very large and forms a quasi-continuum, but a major one in the quantum regime where this 
number is small. 

In recent years, the Gaussian Binning (GB) procedure (see the next section) was fre- 
quently used to include Bohr quantization of product internal motions in QCT calculations, 



removing thereby to some extent the previous drawback 14H21J . However, this procedure is 
numerically inefficient for larger than four-atom systems. This is a severe limitation owing 
to the increasing number of polyatomic processes under scrutiny today. 

Recently, a particular implementation of the exit-channel corrected phase space theory 
of Hamilton and Brumer was proposed, which avoids the use of the GB procedure while 



taking into account the quantization of product vibrations [22J, [23j. In spite of the fact that 
this statistico-dynamical approach represents a significant advance regarding the description 
of indirect polyatomic fragmentations, it is partly based on the micro canonical equilibrium 
assumption of transition state theory and cannot be used for direct processes. 

Here, we propose an alternative equation for the state-resolved absorption cross section 



which is strictly equivalent to the standard expression. From the numerical point of view, 
however, this new equation appears to be much more efficient than the standard one as far 
as inclusion of Bohr quantization in the classical description is concerned. In particular, no 
Gaussian weights are required. Therefore, this formulation opens the way to more realistic 
classical dynamical studies of direct polyatomic photo-dissociations in the quantum regime. 
The alternative equation is derived in section II and its equivalence with the standard 
one is illustrated from an academic example in section III. Section IV concludes. 

II. THEORY 

A. System 

For simplicity's sake, we consider the collinear photo-dissociation of the triatomic 
molecule ABC leading to A and BC. The generalization of the next developments to re- 
alistic systems is straightforward (though obviously heavier). 

The space coordinates of the problem are the usual Jacobi coordinates R, the distance 
between A and the center of mass of BC, and r, the length of BC. The conjugate momenta 
of R and r are P and p, respectively, u = (R,r), p u = (P,p) and T = (R,r,P,p) are the 
configuration, momentum and phase space vectors, respectively. The kinetic energy is given 
by 

P 2 p 2 
T=— + i- (2.1) 

2/i 2m V ; 

where // is the reduced mass of A with respect to BC and m is the reduced mass of BC [5[. 
We call V the potential energy of the photo-excited molecule. Its minimum is supposed 
to be zero in the separated products. E is the energy available to the products. 

B. Standard formulation 

Before the optical excitation, ABC is in the stationary state *£?(u). In the classical 
description of the absorption cross section, the corresponding distribution of the phase 



space states is identified with the Wigner density [5] 



i. 



Mr) 



{nhy 1 [ dv ^*(u + v) *(« - v) e 2ivpu/h . (2.2) 



Moreover, the excitation is assumed to instantaneously occur at time 0. Hence, each quadru- 
plet r determines the initial conditions of a trajectory eventually reaching the product chan- 
nel. pw(T) is thus also the phase space distribution of the photo-excited molecule at time 
0. 

Save for an unimportant normalization constant, the standard expression of the state- 
resolved absorption cross section reads [5j 

a(cu,n) = u J dT p w {T) d(uf 8{H - E) S(x{T) - n). (2.3) 

bj is the frequency of the photon (in unit of 2ir), d(u) is the component of the transition 
dipole function in the direction of the electric field vector, 

H = T + V (2.4) 

is the classical Hamiltonian, n is the vibrational quantum number of BC and x(T) is the 
value of the vibrational action in the separated products. 

In practice, the delta-functions in the previous equations are approximated by functions 
normalized to unity. The delta-function constraining the energy is usually replaced by 
a narrow bin as compared to the width of the absorption spectrum though a Gaussian 
function can equally be used. The delta-function constraining the action was generally 
replaced by a unit-size bin in the past, a procedure called standard binning (SB) or 
histogram method [l2|, [l3j. However, this procedure may lead to unrealistic predictions 
when a small number of quantum states are available to the products, as shown in the 
recent years. As previously stated, a possible alternative to the SB procedure is the GB one 
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21| , where the delta- function constraining the action is replaced by the Gaussian function 

-x 2 /e 2 

G{x\e) = e ——-, (2.5) 



e being usually kept at ~ 0.05. This Gaussian is ~ 10 percent wide, meaning that ~ 90 
percent of the trajectories do not contribute to the dynamics. Consequently, ~ 10 times 
more trajectories have to be run in order to get the same level of convergence of the results 
as with the SB procedure. However, we are dealing here with a very simple system involving 
only one vibrational degree-of-freedom. With N degrees, ~ 10^ times more trajectories 
have to be run. Since in polyatomic processes, N can easily be 10, it is quite clear that the 
amount of trajectories necessary to converge the calculations within the GB procedure is 
just prohibitive. We now consider a simplifying change of variable to go round this problem. 

C. A simplifying change of variable 

In a first step, we assume that BC is rigid. The phase space variables are thus reduced 
to R and P. 

For a repulsive potential energy V exponentially decreasing with R, phase space tra- 
jectories have typically the shape of the path represented in Fig. [TJ This path, defined 
by 

H=^ + V = E, (2.6) 

2/z 

comes from infinity with a constant negative momentum Pi = — (2/zi?) 1//2 , slows down due 
to the repulsive wall, touches a turning point and subsequently follows a symmetric path 
with respect to the i?-axis up to infinity. 

We now consider a large value Ri of R for which V is negligible (see Fig. [1]). Any point 
along the previous path can then be defined by time t such that at t = 0, R = R4. It is clear 
that H and t uniquely define any point of the phase space (R, P). 

The first order developments of dR/dt and dP/dt in terms of t are found to be given by 

dX/dt = Xi + Xit, (2.7) 

X = R,P, with 

Ri = —, (2.8) 



p = - 



dV_ 
dR 



R—Ri 



(2.9) 



Ri 



ldV 
/J, dR 



R—Ri 



and 



d 2 V 



dR 2 



Pi. 

R=Ri LI 



(2-10) 



(2.11) 



Eqs. (|ZSJ) and flZSJ) are Hamilton equations. Eq. f[27T0l) is just (EH]) divided by /i. Eq. 
(12.1 ip is obtained from deriving the right-hand-side of (12. 9 p by R and multiplying by the 
velocity (I2.8p . Moreover, from the first order developments 



X — Xi + Xjt, 



X = R,P, and the fact that 



(2.12) 



P; = -[2,U(#-V)] 



1/2 



(2.13) 



is a function of H, one deduces the two identities 



dR/dH 



and 



P; 



/' 



OP/dH = ^ 



(2.14) 



(2.15) 



From Eqs. ([2TTj) -( l23T]) . fT2~T4j) and (121511 . we quickly find that the Jacobian of the 
transformation (R, P) — > (t, H) is equal to unity, i.e., the transformation is area preserving. 
We thus have 



dRdP = dtdH. 



(2.16) 



As we never used the fact that at Ri, V and its derivatives are zero, the previous result is 
general, i.e., valid anywhere along any trajectory. 

D. Alternative formulation 

Since we are still supposing that BC is frozen, we consider the absorption cross section 
c(uj) given by Eq. (12. 3ft without the last delta-function. From Eq. (12.16J1 . we can replace 
the volume element dT = dRdP by dtdH in Eq. ( 12.31) and perform the trivial integration 
over H leading to 

/+oo 
dtp w d 2 . (2.17) 

-oo 

The time-dependence of pw and d is implicit. 

If we relax the constraint on the rigidity of BC, Ri, t and H are not sufficient to specify 
the dynamical state of ABC. The initial values r$ and Pi of r and p are also necessary. After 
some steps of algebra analogous to the previous ones, we again arrive at the conclusion that 
the transformation (R,r,P,p) — > (ri,Pi,t,H) is unitary, and the state-resolved absorption 
cross section (12.31) turns out to be given by 

a(u, n) = uo / dt dr t dpi pw d 2 5(x — n), (2-18) 

the dependence of pw, d and x on t, r, and pi being implicit. Here, the value of the initial 
momentum P is given by 

Pi = -[2fi(E-E n )]^ 2 (2.19) 

where E n is the vibrational energy in the state n. 

We can also perform the usual change of variable (ri,Pj) — > (q, fix) where q is the angle 



conjugate to x. Since this transformation is canonical 24], the volume element dridpi in 



the above integral can be replaced by hdqdx. Integration over x is then trivial and leads to 



the central result of this note: 

a(u,n) = hu / dtdqpwd 2 . (2.20) 

The numerical efficiency of this expression is obvious: trajectories are now started from 
the products toward the interaction region with x kept at n, i.e., in such a way that Bohr 
quantization conditions are exactly satisfied ; Gaussian statistical weights turn out to be 
unnecessary here, in contrast with the standard approach where the final actions are not 
controlled. 

The statistico-dynamical method for indirect photo-dissociations evoked in the introduc- 
tion shares the same feature. 

Just as a(cu,n) is expressed in terms of inelastic scattering states in the quantum 
treatment [5j , it is expressed in terms of inelastic trajectories in Eq. (j2.20p . This expression 
is thus closer to the quantum one in spirit than Eq. (12.31) which involves only half collision 
paths. 

III. APPLICATION TO A MODEL SYSTEM 

We shall now apply Eqs. ( 12.31) and (12.201) to the well-known Secrest and Johnson 
potential energy [25] given by 

V = -mw 2 r 2 + e a{r ~ R) . (3.21) 

The free BC diatom is thus a harmonic oscillator. For the sake of simplicity, we kept w, 
[j,, m, h w and a at 1. We assumed that the Wigner density was that of two uncoupled 
harmonic oscillators in the lowest vibrational state. Its expression is of the type 

Pw (T) = G(R\a R ) G(P\l/a R ) G(r\a r ) G(p\l/a r ) (3.22) 

where we recall that G is the normalized Gaussian function defined by Eq. ( 12. 51) . a R and 
a r were both kept at 1. Generating T from pw{T) by Monte-Carlo sampling led to total 
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energies H mostly in the range [0, 12]. We kept E at 6, a value for which the density is still 
large and 6 vibrational levels are available. 

We start with the practical evaluation of Eq. (12. 3p . The delta- functions constraining 
the energy and the action were both replaced by the Gaussian functions G(x\e e ) and 
G(x\e a ) with e e = 0.15 and e a = 0.05. T, selected according to pw(X) as previously stated, 
determines the initial conditions of a trajectory calculated using a Runge-Kutta integrator 



26[ up to the separated products defined by R = 10. A batch of iV equal 10 million 
trajectories was run. A time increment r of 0.01 was used for the numerical integration. 
The MC expression of cr(u, n) reads 

1 - 
<r(w, n) = - J2 G ( H * ~ E \^e) G{x{T i ) - n\t a ) (3.23) 

i=l 

where Hi and Ti are respectively the values of H and T for the i th trajectory. 

As far as Eq. (J2.20P is concerned, trajectories were started from a large distance Ri = 10 
with Pi given by Eq. (12.191) and 

E n = n+^. (3.24) 

q was randomly selected in the range [0, 2ir] and r^ and Pi were deduced from the identities 

ri = (2n + l) 1/2 cosq (3.25) 

and 

Pl = _(2n + l) 1/2 sinq. (3.26) 

The analogous (far more complex) transformation for polyatomic processes is given elsewhere 



22]. 



The trajectories were then integrated within the Jacobi coordinates until they recrossed 
the line R = Ri in the product direction. A^ = 200 trajectories were run for each value of 
n. The MC expression of a(u,n) is 

2 N Ni 

i=l fe=l 



The time increment r has been previously defined, Ni is the total number of time steps along 
the i th trajectory and p% is the value of pwiX) for the i th trajectory and k th time step. 

The two sets of predictions are represented in Fig. | Asa matter of fact, they are in 
excellent agreement, illustrating thereby the equivalence between Eqs. (12.31) and (I2.20p . 

With the standard method, 10 million trajectories were found to be necessary in order to 
recover the results obtained with only 1200 trajectories when using the new method. It is 
however clear that the first approach provides the cr(u, n)'s for all the available values of E, 
not only for E = 6. When taking this fact into account, one arrives at the conclusion evoked 
in the introduction, i.e., the amount of calculation is one order of magnitude less with the 
new method than with the standard one, for there is only one vibrational degree-of-freedom. 
In the case of 10 degrees, however, the numerical saving is just amazing. 

IV. CONCLUSION 

A new classical equation for the state-resolved photo-dissociation cross section has been 
presented. This equation is strictly equivalent to the standard one, but is numerically 
much more efficient as far as including Bohr quantization of product internal motions in the 
calculations is concerned. Therefore, this new expression opens the way to more realistic 
classical dynamical studies of direct polyatomic photo-dissociations in the quantum regime 
where only a few states are available to the products. 
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Fig. [T] Typical trajectory in the phase space plane (R, P) for a repulsive potential V. 
The trajectory comes from the separated fragments with a negative momentum 
Pi, slows down due to the repulsive wall, touches a turning point (red dot) and 
follows a symmetric path in the upper plane. An arrow indicates the direction 
of motion. Time t = corresponds to passage at R = R{, a value for which V is 
vanishingly small. 

Fig. [2J State-resolved absorption cross section in terms of the vibrational quantum num- 
ber n of BC. The red and blue curves are found from Eqs. (12.31) and (j2.20p 
respectively. 
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